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i" - !" Abstract 

cu. 

This paper deals with a nonlinear system of partial differential equations modeling 
a simplified tumor-induced angiogenesis taking into account only the interplay between 
tumor angiogenic factors and endothelial cells. Considered model assumes a nonlinear flux 
at the tumor boundary and a nonlinear chemotactic response. It is proved that the choice 
of some key parameters influences the long-time behaviour of the system. More precisely, 
we show the convergence of solutions to different semi-trivial stationary states for different 
range of parameters. 
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1 Introduction 

Angiogenesis is a physiological process involving the new vessels sprout from a pre- 
existing vasculature in response to a chemical stimuli. Angiogenesis is an important ingre- 
dient of a processes like development, growth and wound healing. However, angiogenesis 
is also induced by tumoral cells. In this paper we consider a model of tumor-induced an- 
giogenesis that was proposed in [5J. Actually, in the above mentioned model some factors 
influencing angiogenesis are neglected to keep the model simple but sufficiently interesting 
from the analytical point of view. We refer the reader to [T2] as a source of information 
about the progress in mathematical modelling and biological knowledge of angiogenesis 
process. We focus our attention on two key variables: the endothelial cells (ECs), denoted 
by u, and the tumor angiogenic factors (TAF), denoted by v. We assume that (ECs) 
that form the blood vessels wall are induced by the (TAF), factors that are generated 
by the tumor, to migrate chemotactically towards the tumor. We assume that the (ECs) 
and the (TAF) fill in a bounded and connected domain C H d with a regular bound- 
ary 80. In particular, neither the existence of extracellular matrix nor the activity of 
metalloproteinases is considered. But, what was new there, nonlinear flux of TAF on the 
tumor boundary was taken into account. The reason was that since ECs are supposed to 
react chemotactically to the TAFs, generating the large gradient of TAFs on the boundary 
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would probably make the tumour more dangerous. The aim of [5] was to study the inter- 
play between the density of ECs and TAFs dependently on a parameter \i measuring the 
strength of the flux on the tumor boundary and the nonlinearity V measuring nonlinear 
response of ECs. In |2] the qualitative features of a model were studied in a local sense. 
We mean by that the local stability of steady states which were proven to exist in [5] . We 
complete the studies taken in [S] by analyzing the global stability of steady states. We 
shall prove the asymptotic convergence of solutions for different values of /x. To be more 
precise, we consider the case 

d0 = r 1 ur 2 , 

where ri n T 2 = and Tj are closed and open sets in the relative topology of 80. We 
suppose that r 2 is the tumor boundary and Ti is the blood vessel boundary. Our parabolic 
problem reads. 



Au 



vt — Av = 
du dv 

dn 
du 



dn 
0, 



-div(V{u)Vv) + lu 
-v — CUV 





dv 
dn 



dn dn 1 + v 

u(x, 0) = uq(x), v(x, 0) = vq{x) 



in x (0,T), 
in x (0,T), 

on Ti x (0,T), 

on T 2 x (0,T), 
in 0, 



(1) 



where < T < +oo, /, u € H, c > 0, 



V 6 C X (R), V > in (0, oo) with V(0) = 0; 



and no and vo are given non-negative and non-trivial functions. In Theorem 3.1, 
Theorem 3.8] the existence and uniqueness of global-in-time bounded regular solutions, 
provided initial data are nonnegative and V G L°°(0,+oo) is shown. Moreover in [U 
Section 4] the existence of two semi-trivial steady-states (A,0), A > and (0,0u) is shown 



provided u > u\ (see also [14]) 
eigenvalue problem 



where ui is the principal eigenvalue of the boundary 



— Av + v 

^ = 
dn 

dv 

dn- = flV 



in Q, 
on Ti, 



on r 



2- 



Furthermore, results concerning the linearized stability around the semi-trivial solutions 
to ([T]) are proven in [5]. 



First models of tumor induced angiogenesis that we are aware of are considered in [3] 
(see also [10] for a more elaborated model). A reduced model proposed in [10] is studied 
in [7]. The local stability of the homogeneous steady-states in one dimensional domains 
is shown there. In all the mentioned papers the boundary conditions are either zero Neu- 
mann or no-flux. In [6] the stationary problem of (HJ with linear flux for v is studied. 
Finally let us mention [9] where the authors study the local solvability of a system of 
partial differential equations with a nonlinear boundary condition and a chemotaxis term. 
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The aim of this paper is to analyze the global stability for positive initial data. In 
particular we show global stability for some range of parameters (A, /x) for which even the 
local stability is not known. 

It should be pointed out that the results of this paper could be extended even to more 
general forms of V as soon as 

IK*) Hoc < C, for t > 0. (3) 
Observe that, if the above inequality holds, then the parabolic regularity asserts 

IKi)IU<C 

for any t > and by [SJ Theorem 3.1] the solution is global and regular. In particular, 
when V is bounded in the L°° norm (see [5]) then ([3]) is satisfied. 

2 Preliminaries 

For the reader's convenience we collect here some results of interpolation theory and 
its applications to parabolic problems that will be used throughout the paper. 

a) Let Eq, Ei two normed spaces, we can define the real interpolation functor, denoted 
by 

(E ,Ei)e, p , O<0<1, l<p<+oo, 

(see for instance |13t Def. 22.1]). During the paper we will use the following property 
of the real interpolation functor (see |13|, Lemma 25.2]): 
If (Eq, Ei)$ iP is a Banach space then 

3C > such that ||a;||( BoiBl ) e>j) < CH^H^""^!!^ Vx E E n E t . 

In the context of fractional Sobolev spaces this inequality reads, see [U Theorem 7.2] 

\\u(t)\\ W m, P <c\\ u (t)\\ e wk , p \\u(t)\\l- e (4) 

for m<k9,9 G (0,1). 

b) Let us consider a parabolic problem with a non-homogeneous boundary condition 

z t + Az = f(t) in n x (0,T), 
< Bz = g(t) on9Ox(0,r), (5) 

z(x, 0) = zo(x), in £1. 

where 

on 

and 

Az := -Az + z. 

We define the space of functions 



{z G W S ' P (Q.) : Bz = 0} ifl + l/p<s<2, 
W S > P (Q) if -1 + 1/p < s < 1 + l/p, 

(W~ s ' p '(n))' \i -2 + l/p < s <-l + l/p. 
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It is known that (A,B), as being in separated divergence form (see [U pg. 21]), is 
normally elliptic. We denote by A a _\ the Wg° 2 ' p -realization of (.4, B) (see [U pg. 
39] for the precise definition). Since (A,B) is normally elliptic then A a _\ generates 
an analytic semigroup [TJ Theorem 8.5]. Moreover, if 

(f,g) G C((0,T);Wl a - 2 > p (n) x wf l ~ 1/KP (Xi)) 

for some T > and 2a £ (1/p, 1 + 1/p) then for any t < T we rewrite © by the 
generalized variation of constants formula 



z(t) 



-tA 



'z + 



-(t-T)A a 



^{f{r) + A a ^Blg{r))dT, 



where B c a is the continuous extension of {B\ Ker ^)~ 1 to W 2a ~ 1 ~ 1 / p ' p (dQ,). Since 
[0, +oo) C p(— A*-i) (p is the resolvent set) then by [H Remark 8.6 c)] there exists 
a constant C > 1 such that 



14' 



iot,p ^ C||y4.(j — 1^| 



2a — 2,p . 



(6) 



c) Let a, b, c € L°°(0), the eigenvalue problem 

-Az + a(x)z = Xz in fl, 
on Ti, 



— + b(x)z = 
an 



<9-£ 

— + c(x)z = 
on 



on r 2 . 



has a unique principal eigenvalue (i.e. an eigenvalue whose associated eigenfunction 
can be chosen positive in £7) and it will be denoted by 



Ai(-A + a;M + b;M + c) 



3 Convergence to the semi-trivial solution (/,0) 

In the present section we deal with the convergence to the semi-trivial steady-state 
(1,0). Throughout this section we assume ([3]). A sufficient condition guaranteeing ([3]) is 
the boundedness of V (see [5]). We will use the generalized variation of constants formula 
to estimate v, which is stated in the next lemma. 

Lemma 3.1. Let 7 G (l,+oo) and (5 £ (1,2a). Then, for every r £ (0,t) there exists a 
constant 5 E (0, Re a(A a -i)) (a denotes the spectrum) and 9 = 6(j3) € (0, 1) such that 

\\e-^ A "-iz\\ w ^ < C(t-r)- e e- s ^ \\z\\ w 2c-^ 



for every z € W s a ' 7 . 

Proof. By the choice of (3 we have W[ 
Theorem 7.2] we get 



= W^'T(f2). As a consequence if we apply [H 



e -(t-r)A a -, z \\ wf)i < C\\e-^ A ^z\\ w 2^\\ e \\e-^ A ^z 
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for some 8 G (0, 1). Next we apply (|6]) to the first norm on the right hand side and 
Theorem 1.3.4] to deduce 

where 5 G (0, -Re a(A a ^i)). m 

Lemma 3.2. Let 7 G (1, +00), /3 € (1, 1 + I/7), // e [0) ^1) ari ^ < 5 < p < a(/x) u>/iere 
a(/u) is defined as 

a(/i) := Ai(-A + l;J\f,J\T- /x). 
Then, there exists C > suc/i i/iai, /or i > 0, £/ie v-solution to (OP satisfies 

v(x,t) <Ce- pt V(x,t) G H x (0,+oo), 

||«(*)||^,7 < C(l + t- 9 )e- ft ||«o|| 7 , 

w/jere 6> = 6>(/3) G (0,1). 

Proof. A solution to the problem 

w t -AH! + w = in Q x (0,T maa! ), 

■^- = on Ti x (0,T mm ), 
dw 

-g^=IJLW on T 2 x (0,T ma3; ), 

w(x, 0) = vo(a;) in f2, 



(7) 



is a supersolution to the -u-equation of (HJ, therefore v(x,t) < w(x,t). Since, for suffi- 
ciently large M, W = Me~ pt <pi, with epi a positive eigenfunction associated to a(/i), is a 
supersolution to ([7]), the pointwise estimate in the claim of the lemma follows. For the 
second one we pick 

f(t) := -cu(t)v(t), 

{0 on Ti, 

v(t) r 
ix tt on 1 2 • 

Taking the -norm in a generalized variation of constants formula for v and using 
Lemma 13. II we obtain 

-tA a -i-. 11 „ , f \\„-{t-T)A a 



\W)\\we«< l|e~ tAQ -^o|lH/^ + / We-^^UiT) + A Q ^B c a g(r))\\ w ^ 

t 

< C (e- st t- e \\v \\ w 2 a -2,y + (t - T y e e- 5 ^\\f(T) + A a ^B c a g(T)\\ w 2 a -2,-,d^j . 
Next, we estimate the last term in the above inequality using the fact that 

A a ^B C a G £(^2a-l-l/ 7 ,T '(0^ w2a-2o , ^ 

and the continuous embeddings 
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Therefore, we get 

Ht)\\ w ^<e- 5t t- e \\v \\ y + Ce- 5t re 5T (t-r)- e (||/(T)|| i7(n) + || 5 (T)|| i7(an) )dr. (8) 

J o 

Observe that by ([3]) and the first part of the Lemma we have 

\\f{r)\\L^ m <C\\v\\ Loom <Ce~P\ 

h(T)\\ L -i{dCi) < \\v\\l<x>(8CI) < Ce~ pr . 
In view of the above bounds, (|SJ) yields 



IK*)||w^ < Ce~ st t~ e \\v \\j + Ce~ 5 *|[u || 7 / e (5 - p)T (t - r)- e dr. 

■/ o 

/■OO 

Next, by the choice of 5 and p, / e^ s ~ p ' T (t — t)~ 9 cIt = C < +oo and the Lemma follows. 

J o 

■ 

Our purpose is to show that u converges to steady states. To this end we treat sepa- 
rately the cases A = 0, A > 0. 

3.1 Case A = 0. 

Lemma 3.3. Let t > and y G C x (t, +oo) fl L'(t, +oo), y' G L 1 (r, +oo). Then 
lim |y(t)| = 0. 

t— >-+oo 

Proof. By the assumptions of the lemma we observe that for any k > 

lim / m 0( S )| + |y'( S ) |)ds = 0. (9) 



Let us assume that lim \y(t)\ ^ 0, then there exists a sequence {t n } n£ -^, t n — > +oo, such 

f— !-+oo 

that 

\y(tn)\ >C>0, 

for all n > hq. We pick 9 G (0, k], then for any e > 

rtn+6 rtn+k 

\\y(t n + 0)\ - \y(t n )\\ < \y(t n + 9) - y(t n )\ < \y'(s)\ds< \y'(s)\ds<e 

by ©. Therefore \y(s)\ > C/2 for all s G [t n , t n +k], n > no- The last assertion contradicts 
the fact that 



rtn+k 

lim / \y(s)\ds = 0. 



In the following lemmata (ii,i>) is a solution to (pQ). 
Lemma 3.4. Let A = and t > r > 0, then it holds 



r 2 1 + u Jn Jn Jn 



A* "7^7 + u = «W- «(t). (10) 
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Proof. Integrating the u-equation of ([I]) yields 

Ut f V(u) — 

Jdn \dn dn t 

So, integrating the last expression in time between r and t we get the result. 
Remark 3.5. By Lemma \3.J\ we see that for any t > r 



u < \\u(t) 



II- 



n 



Theorem 3.6. Assume that < fx < fi\ and A = 0, then 

lim \\u(t) \\w m * = 0) 

t— s-+oo 

for any m < 1 anc? p > 2. 

Proof. On multiplying the u-equation of ([!]) by u and integrating in space we obtain 



V(u)uv 



& 1/ = X (H V "' 2 + V{U) V " ■ V " - " 3) - " /r a T + . 
<(«-!) /|V»|' + C7 W /|V.|»-m/ 

Therefore, we infer 

A ( u 2 + (l-e) [ \Vu\ 2 <C(e)\\v\\ 2 wl , 2 , 

and after integrating in time, thanks to Lemma 13.21 we arrive at 

u(t) 2 - [ U {rf + {l-e) f [ |Vu| 2 <C(e) f\l + s' e ) 2 e~ 2Ss \ 
i Jn Jt Jn Jt 



u 



vol 



In particular we deduce that for t > r 



|Vd 2 < C. 



By Lemma 3.8] we find a bound < C, therefore 



_d_ 
2dt 



[u 2 <C [ \Vu\ 2 + C{e)\\v\\ 2 wi , 2 +cJ' Y^l + cf 
Jn Jn Jr 2 1 + v Jn 



u 



Thanks to (TO . for t > r 



2dt 



< C. 



Finally, Remark 13.51 and (|12p together with Lemma 13.31 entail 



(11) 



(12) 



lim Ki)|| 2 = 0. 
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Also thanks to ||w(t)|| c ^ < C for all t > we obtain 



lim \\u(t) |L = 
t->+oo ^ 



for any p > 2. Next we recall that by O Lemma 3.7] for any 2/3 € (k, 1) we find a bound 
on the Xp norm of u, where Xp is a usual fractional space connected to a semigroup 
approach to parabolic equations, see [8]. Next, due to the fact that 2/3 € (k, 1), we infer 
from the embedding Xp (see for instance [SI Theorem 1.6.1]) that for all k < 1 

and p > 2 

\\u(t)\\ W k, P < C. 

Next, (gD entails 

Ki)||w"».p < C||n(t)||^ fc , p ||n(t)||i- e . 

Therefore, it holds 

lim \\u(t)\\ W m, P < C lim \\u(t) = 0. (13) 



Remark 3.7. Lei pomt out t/iai if we pick m such that m — d/p > then W m ' p (fi) is 
embedded in C(Q). 

3.2 Case A > 0. 

Assume that there exists 5q and to such that 

u(t) > 6 > (14) 

for t > to > 0. Next, we examine the long time behavior for u under the hypothesis ()14p . 
In the sequel we shall give sufficient conditions on V(v) implying (I14p . 

Theorem 3.8. Let < fj, < m and assume the the hypothesis JHJ)) is satisfied, then there 
exists 9 > suc/i that 

\\u(t) - \\\ w ™,p < Ce~ et , (15) 
for all t > to and any m < 1, p > 2. 
Proof. On multiplying the ^-equation by u — A we have 



d 

2dl 



j\ u -\f= - [ |Vu| 2 + / V(u)Vv • Vu - n I yp&fa-X)- f 
1 r in/ii2 /• 



u(« - A) 2 



< -I [ |Vul 2 + ^i / |vd 2 



2 Jn 2 j n 

J T (Tt^f) - J n u(u ~ x)2 - 

(16) 

Having in mind that (1 + t>) 2 > 1, the hypothesis f)14[) and the Sobolev trace embedding 
we get 

^ f (u - A) 2 + 25 / (« - A) 2 < (7^11^,2 + mC|M| w i,2. (17) 
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By Lemma 13.21 we can deduce 

IK*)- A||! < Ce~ 6lt 

for < Q\ < min{2<5o, /3}. At this point we can argue exactly as in the end of the proof of 
Theorem 13,61 Namely, by the bound on u in L°° we infer the bound on the IP norm of u, 
p > 2. Next, we use the estimate of u in W ,p , k < l,p > 2, coming from [5., Lemma 3.7], 
in order to conclude (fT5l). ■ 



In the rest of this section we give sufficient conditions on V implying (|14|) . Actually, 
only the behavior of V around zero matters. Roughly speaking we require a superlinear 
growth of V in the neighbourhood of zero. From now on we assume that there exist 
C, 5 > 0, k > 1 + d/2, j > d/2 such that 

< V(s) < Cs ko , \V'(s)\ < Cs j (18) 

for all s G (0,5). 

Remark 3.9. The condition (]18p is satisfied, for example, for functions 

V{u) = 

with a > 1 + d/2. 

Next we introduce some notation that will be of importance in the proof of (I14[) . 
Moreover we formulate a lemma which we need in the main part of the proof of (|14p . Let 
f(S),g(6) be defined in a following way: 

f(6) := sup V 2 (s), 
se(o,<5) 

g(S):= sup (2(s-<5)^y / (s) 2 + 2y 2 (s)). 
se(o,<5) 

Lemma 3.10. Assume that 173)) holds. Moreover, for some D,fj, > 0, r/ > 1, e and C(e) 
are given by 



2nD' K ' 26 2r >' 

Then, if 5 > is small enough, the following conditons are satisfied simultaneously 

V 2 (s) 

C(e)—^S<X-6 (19) 
s 

for s G (0,(5) ; 

C(%(5) < 1/2 (20) 

and 

fW < — . (21) 
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Proof. Thanks to (I18p . we have 

f(8)D = sup V 2 {s)D < C5 2ko D. 

se(o,s) 

Hence, for n < ko and 8 sufficiently small ([2T]) is satisfied. Next, owing to (fT8|) . we observe 
that 

s 

Thus, ()19p can be assured for rj < ko and 8 small enough. Moreover, it is straightforward 
to see that ([20]) is also satisfied for 1 < r\ < min{/co, 1 + j}- ■ 

Lemma 3.11. Assume that < fi < fj,\ and that (|18p is satisfied then (|14p holds. 

Proof. Let (5 > be a fixed constant defined in (|18|) . Given a function /, we define the 
negative part of / as a nonpositive function as follows 

/_ := min{/,0}. 

Our purpose is to show that ||(it — <5) — (t)||oo < 8/2 for every t > to which implies (JHJ). 
In order to obtain the previous estimate we multiply the u-equation by (u — 8)- and we 
integrate in space to obtain 

+L(f^«S0 ( "-w„" (A -"> ( "-^ 

= - f |V(u-<5)_| 2 + / V(u)Vu- V(u-<5)_ 

WitW— — («- S)- + / u(A-m)(u-£)_ 
1 + v 



JT 2 1 "f " 

/ |V(u-<5)_| 2 + / f(ti)Vt)' V(u-<5)_ 
Jn Jn* 



/// V(«)(«-(5)_ + / u(A -«)(«-<*)_, 



r 5 1 +f ./<> 



where 



Consequently, 



fi<5 := {x G ^ : < 5} , := {x G 1^ : < 5}. 



± [{u-8f_< ( e -l) f \V(u-8U 2 + C(e) [ V 2 (u)\Vv\ 2 

-li [ —^V(u)(u-8)-+ [ u(X-u)(u-8). 

< (e-1) / jV(u-5)_| 2 + C(e) sup T/ 2 (s) / |V 

s6(o,<5) Jn 

-u / — —V(u)(u-6)- + / «(A -«)(«- <J)- 
Jr* 1 + u Jn 



v\ 2 
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Previous inequality can be rewritten in terms of f(5) denned before Lemma 13.101 as 
"' r {u-8f_< (e- 1) f \V(u-S)^\ 2 + C{e)f{S) [ |V«| 2 + /i? 1 



+/iC(e) / V(v)\u-S)t+ / u(A -«)(«- (5)_. 
■/r 2 ./n 

Thanks to the Sobolev trace embedding W ' 2 (0) ^4 L 2 (<9$7) and having in mind that 
(v + l) 2 > 1, we arrive at 

V{u) 2 {u-S) 2 _< C ( [ V 2 {u){u-5) 2 _ + 
r 2 \Jn 

<> 

2 



Therefore, we obtain 
d 

2dt 



+ / (2(n-5) 2 _F'(n) 2 + 2y 2 (u)) |V(u - <5) 



[ (u-S) 2 _<(e-l) [ \V(u-S)^\ 2 + C(e)f(S) [ \Vv\ 2 + Ce\\v\\ 2 m , 
Jn Jn Jn 

+C(e)( f V 2 {u){u-S) 2 _ + I (2(u- S) 2 _V {uf + 2V 2 {u)) \V{u - 5)- 
\Jn Jn 



+ / u(A -«)(«-*)_. (22) 
Jn 

In view of the nonnegativity of u we have 

-6<(u-S)-. (23) 
Owing to ()23p . from (|22[) we see that (g(<5) was defined before Lemma l3.10p 

^7 / (u-«£ <(e+C(e)0(J)-l) / |V( U -5)_| 2 

+ (C(e)f(S) + /ie)||«||^ lia + j^u{u - f A - u - C(e)^-^ 
Due to the nonpositivity of (u — 5)- and (I20p we have 



(24) 



u(u-<5)_ (A-u-Cfe)^^ <0. (25) 
n V u J 

By the Hopf lemma and zero Neumann data on the boundary for u we see that there exists 
S\ such that u(to) > Si. Hence choosing 5 < Si and using (fT9l) . ([25]) and Lemma 1331 we 
infer from (|24|) 

||(« - ^-(Olll < (2C(e)/( < 5) + 2^?)C(/3), 
for i > to > 0. We shall show that 

\\(u - SUt)\\ 2 < S 2 \ (26) 
11 
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for some 77 > 1. To this end notice that choosing e = C(e) = 1/2, we are in a position to 
apply Lemma [3.101 with D = C(/3). As a consequence, for e as it is chosen in Lemma [3.101 
PI), CE]), (ESQ) and 

5 2r > 

2»eC(/3) < — 

are satisfied simultaneously. Hence (I26p is shown. 

Next we use interpolation between L p spaces, (|26p and (|23p to obtain 



\\{u-5U\ 2/ei <||(u-*)_||f||(«-«)_||^ 

< §010! §1-9! _ jl+(a-l)0i_ 

Applying (|1J), we infer 

|| (« - 5)_ 11^,2/^ < C||(u - SU\ e wh2/ei \\(u - SU\ljl ^ C ^ u ~ V-Wlfr, 

the last inequality being a consequence of the uniform bound of L°° norm, see [U Theorem 
Lemma 3.8], and [H Theorem 15.5]. Picking up 0\ such that 

°- d -Y>° ( 27 ) 

we make sure that W e ' 2 l dl (Q) <^ L°°(0). Consequently, 

\\(u- 5)^ <C 2 5^ 1 - e ^ a - 1 ^. 
Next, we notice that choosing a > 1 + | we make sure that 

l<(l-^) (l + («-l)0i). 

Hence, choosing 9 close enough to we see that (1 — 6>)(1 + (a — l)#i) > 1 and upon 
taking 5 small enough we obtain 

\\(U - SUt)\\oo < ^, 

for t > to > 0. The Lemma is proved. ■ 

4 Convergence to the semi-trivial solution (0, 6> M ) 

Through this Section we additionally assume that there exist constants < c m < Cm 
and a > 1 such that 

c m s a < V(s) < C A fS Q for all s € [0, |M|oo]. (28) 

Remark 4.1. Let us observe that when V'(0) ^ and (0j holds, then I128\) is true for 
a = 1. Moreover if V e C k for k > 1 with V k (0) 7^ and F>(0) = /or j < fe, ^en fig)) 
holds true for a = k. 

In the following Theorem, we eliminate the restriction on fj, of Theorem 13.61 However, 
we require the additional condition (f28|) on V. 
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Theorem 4.2. Let A = and assume (28i). then 



lim ||u(t)||vK m .p = 0, 

t— ¥ + CO 



for any m < 1 and p > 2. 



Proof. On the one hand, we multiply the it-equation of ([I]) by u and we integrate in the 
space variable to obtain 

2dtJ n * "L l + v (29) 

V [u)uv y ' 



[ (-\V u \ 2 + Vv-V<p(u)-u 3 )-» I 
Jn ' Jr 2 1 + v 



with 

v (u) = / V(s)ds. 
Jo 

On the other hand, we multiply the ^-equation of (pQ) by (p(u). Integrating in space, we 
obtain 

/ Vv-Vip(u) = — [ ip(u)vt + fJL [ *f — [ vip(u) — / CUV(f(u). 
Jn Jn Jv 2 l + " in Jn 

Inserting the above equality into (|29p we have 

d f 2 f / i. - i9 . / \ / \ ^\ f v 



2dt 



(— iVul + ip(u)(— vt — v — cuv) — u ) + a / ((p(u) — V(u)u). 

Jr 2 1 + V 

(30) 

Next we estimate vt- Multiplying the v-equation by vt and integrating over Q we see that 
where 

0(v) := [ V -^—ds. 
K ' Jo 1 + s 

Therefore, by the uniform bound of v in C(Q) we deduce 

~ + f |VH 2 + — [ 8(v)<M f u\ 

4 J n 1 2dt J n 1 1 2dt J n dt J T2 v ' ~ J n 

After integrating over the interval (r, t) we find, by Lemma |3.4| that for t > r 

f [ «? < a 

Jt Jn 



(31) 



Next, by {28} we obtain from (J30J) that 

^ f - 2 ^ /" IV7 |2 i /" ^\2, /" 2, ../" uC M U a ~ l " 



2dty n y n y n 7 Q y r2 (i + u)(a + i) 



<_ / | V n| 2 + max V 2 ( S ) / u 2 + / W 2 + /" 



(32) 



" + 1 Jr 2 1 + « 
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By Lemma 13.41 and (j28j) we get 



vu 



< C 



(33) 

for t > r. According to (|33|) and (|3ip we find upon integration of (|32p over the time 
interval (r, t) that for t > r 

/•* r 

|Vd 2 < C. 



From the last estimate, a similar argument to the one used previously yields 



d 

— I w 

dt 



< C 



for t > t. Thus, by Lemma 
Finally, we can infer the result arguing as in the end of the proof of Theorem [3? 



hm Ki)|| 2 = 0. 

t— >+oo 



Next we prove a lemma which we will use in the proof of Theorem 14.41 As a by-product 
of the following lemma we learn a qualitative information that v is bounded away from 
for times large enough. We shall obtain a lower bound on v by considering a subsolution 
to an elliptic problem which is also a subsolution to a second equation in ([!]). 



Lemma 4.3. Let A = and \x > If the condition I128\) is satisfied then there exist 
constants c\ , r > such that for t > r 



v(t) > C\. 



(34) 



Proof. Let k S (fii,[i). Since Ai(— A + 1;N;N + b(x)) is increasing with respect to b 
(see [U Proposition 3.3]), we have 

Ai(-A + l-N-N — y) < Ai(-A + 1; AT; AT - k) < A x (-A + 1; M; M - W ) = 0. 

Therefore, there exists e > such that 

Ai(-A + 1;JV;JV - k) = -ce i.e. Ai(-A + 1 + ce ;M;M - k) = 0. 

Let tpi be the positive eigenfunction with ||v?i||oo = 1 associated to the above eigenvalue 
i.e. if \ satisfies 

-Atfi + (1 + e c)(pi =0 in Q, 
dipi 



0, 



n 



dn 





ktpi 



on Ti, 

on r 2 . 



By Theorem 14.21 there exists to > such that < u(t) < eo for all t > to > 0. We claim 
that there exists 5 > such that w = 6<pi is a subsolution to 



w t - 
dw 

1- = ° 

m 

yw 



Aw + (1 + cu)w = in x (t , +oo), 
on T\ x (to, +oo), 
on T 2 x (t , +00). 
in Q,. 



w 



dn ^ 1 + w 



w(x,t ) = v(x,t ) 
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Therefore v(x,t) > 5ipi > c\. It remains to prove the claim. By the strong maximum 
principle v(x,to) > c > 0. Thus there exists 5 > such that 5(p± < v(x,to). Moreover, 
choosing 5 > such that fc(l + <5) < /i we make sure that 

dw ^ w 
dn ~ ^ 1 + w 

on T2 x (to, +00). Hence the claim is shown and the lemma follows. ■ 

Now we are in a position to prove the main result of this section. To this end we 
make use of the theorem by Amann and Lopez-Gomez, see [2], stating the equivalence 
between positivity of principal eigenvalue and existence of stricly positive supersolution 
of some elliptic problems (the previous version of this theorem for the Dirichlet problem 
was shown in 

Theorem 4.4. Let A = and assume H281), then 



lim \\v{t) -Buh = °- 
Proof. Let z(t) = v(t) — B^. Then z solves the following parabolic problem 



dn 
oz 



Az - z 



CUV 



On ^(l + t;)(l+0 M ) 
k z(x,0) = v {x) - 8^ 



in x (0,T), 
on r x x (0,T), 

on T 2 x (0,T), 
in 0. 



(35) 



We multiply (|35j) by z to obtain 



d 

2dl 



\Vz\ 2 + n 



cuvz. 



,p 2 (1+«)(1 + ^) 

In order to estimate the right-hand side of ([36]) for t > to, we pick 7 > 1 such that 



7 



< 1 



where c\ is given in ([3 



1 + ci 

For each t > to we consider the eigenvalue problem 



(36) 



(37) 



—Aw + w 

^ = 

dn 

ow 



Xw 



dn (i+v(t))(i + e IJ 



Next, we see that Bn is a strict supersolution of 



in 0, 
on Ti, 

on r 2 . 



(38) 



— Aw + w = 
dw 







M 
)w 



fiJW 



dn (1 + v(t))(l + B p 



in 0, 
on Ti, 

on r 2 . 
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Indeed, 

-A6> M + 6> M = in Q, 

— ^ = on Ti, 

on 

Finally by the choice of 7 (see (I37p ) and Lemma 14.31 we have 

dn + ^ > (l + «(t))(l + ^) ° n 2 " 

Therefore, by [21 Theorem 2.4] we get Ai > 0, the principal eigenvalue of (f38|) . Next, the 
variational characterization of the principal eigenvalue entails 



Jo, Jn Jr 2 \ L 



r 2 



Ai = inf 

ipem(n) 1 2 

n 



Thus, for all tp € H 1 ^) we have 

/V < 7 _1 I |V(^| 2 + 7 _1 /" v? 2 " H / 



(l + u(t))(l + M ) 

In particular, we can apply it in (|36p to obtain the following inequality 



cuvz 



Therefore, there exists M > such that 

4// + (1 -M/J v W/) £M //- 

Integrating the above estimate on the time interval (r, i) we obtain for t > r, 

|Vz| 2 + [ z 2 <C. (39) 



In view of (13911 one infers 

< C 



I - I 

Jr dtj n 



for t > t. Finally, the result follows by Lemma 
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